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ABSTRACT 

We  consider  the  inviscid  stability  of  the  Batchelor  (1964)  vortex  in  a  compressible  flow. 
The  problem  is  tackled  numerically  and  also  asymptotically,  in  the  limit  of  large  (azimuthal 
and  str('amwise)  wavenumbers,  together  with  large  Mach  numbers.  The  nature  of  the  solution 
pass('s  through  different  regimes  as  the  Mach  number  increases,  relative  to  the  wavenum¬ 
bers.  At  very  high  wavenumbers  and  Mach  numbers,  the  mode  which  is  present  in  the 
incompressible  case  ceases  to  be  unstable,  whikst  a  new  “centre  mode”  forms,  whose  stabil¬ 
ity  characteristics  are  determined  primarily  by  conditions  close  to  the  vortex  axis.  We  find 
that  generally  the  flow  becomes  less  unstable  as  the  Mach  number  increases,  and  that  the 
regime  of  instability  appears  generally  confined  to  disturbances  in  a  direction  counter  to  the 
direction  of  the  rotation  of  the  swirl  of  the  vortex. 

Throughout  the  paper  comparison  is  made  between  our  numerical  results  and  results 
obtained  from  the  various  a.symptotic  theories. 
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1.  Introduction 


Ill  i('C('iit  yi'ais  then'  has  Ixh'ii  a  !i;()o(l  (h'al  of  iiitc'n'st  in  tlu'  stability  of  iiicoiiiprcssililc 
swiiliiiu;  vorti'x  type  Hows.  Two  important  apiilications  to  this  an'a  of  rt'st'arch  an-  the 
lircakdown  of  trailing;  liiK'  vortici's  Ix'hind  airi'raft  and  to  tornadoes;  this  class  of  How  may 
also  1)('  applicable'  to  Hows  inside'  turbine's  and  cemiprcssors.  to  which  the'  pre'se'iit  weirk 
weiulel  l)e'  partie'ularly  re'h'vant. 

The'  e'arlie'st  works  in  the'  are'a  of  the'  staliility  of  swirling  vorte'x  Hows  ine'luele'  theise' 
e)f  Le'sse'ii  i.k:  Paille't  (1974)  anel  Le'sse'ii.  Singh  ik:  Pailh't  (1974).  In  the’  fonner  papi'r  the’ 
staleility  eif  the'  Bate-he'Ior  (19G4)  vente'X  was  e'onsieh're'd.  at  Hnite'  Re-ynolds  niimbe’rs  up  tei 
150.  In  the'  se'e-onel  pajie'r,  the'  invisciel  stability  of  this  vorte'x  was  stiulicel  anel  re’ve'ah'el  an 
incre'ase'  in  growth  rate  at  incre'asingly  large  wave'numbers  (with  elisturbane-e's  be'ing  meist 
elange'rems  e-eMinte’r  to  the  dire’e-tiem  of  the’  swirl). 

Duck  k:  Foste'r  119S0)  sheiwe'd  that  for  a  given  wave'immbe'r  a  multijilicity  of  moele's 
e'xists.  Le'ibovie'h  ik;  Ste’wartson  (1983)  and  Duck  (19SC)  consiele're'd  the'  limit  eif  large’ 
wave'iiumbe'r  for  this  preible'in.  anel  sheiwe’el  that  a  Hnite’  (maximum)  greiwth  rate’  was  at- 
taiiK'el.  The’  afeire'ine’iitieme'el  stuelie's  sugge'ste'el  an  uppe'r  anel  leiwe’r  neutral  value’  eif  axial 
wave’numtie'r.  The’  uppe’r  lu'utral  peiint  feir  large  azimuthal  wave’immbe’r  was  tre’ate’el  by 
Ste’wartsein  «k:  Cape’ll  (1985),  wlm  slmwe’el  that  the’  "ring  meiele'"  strue’ture’  eif  the  unstable’ 
meiele’s  pe’rsiste’el  lu’ar  the  upiier  ne’utral  peiints.  Ste'wartsem  Brown  (1985)  ceinsiele'ie’el 
the’se’  uppe’r  ne’utral  jieiints  feir  eirele’r  erne’  azimuthal  wave’immbe’rs,  anel  fenmel  that  in  this 
e-ase’  the’  nmele’s  we’re’  eif  ce’iitre’  ineiele’  type’,  similar  te)  tlieise  femnel  in  a  re’late’el  stuely  em 
swirling  Peiise-uille’  Heiw  ( Ste’wartsein  (C  Breewn  1984).  The’  be’havienir  eif  the  unstable’  meiele’s. 
edeise’  te)  the’  leiwe-r  ne’utral  peiint,  at  large  azimuthal  wave’iiumlie’r,  was  inve’stigate’el  ley  Ste’w- 
artsein  eC  Le’ibeevie-h  (1987),  anel  ele’te’rniine’el  that  in  this  e’ase  the  instability  elisturbane-e’s 
we’ie’  ce’iitre’el  near  the’  axis  eef  the-  veerte-x. 

Me)re’  re’ea’iitly,  visexeus  re’sults  have’  be’e’ii  pre’se’iite’el  at  finite  (but  large)  Re’yneelels 
numbe’rs  by  Kheerrami  rA  al.  (1989)  anel  Kheerrami  (1991).  In  this  latte’r  jeaper,  it  was 
shown  that  aelelitieenal  unstalele'  meiele’s  e’xist,  in  which  visexesity  plays  a  elestaleilising  ixele’. 
The  ■se’  meeele’s  we’ix’  analyse’el  l)y  Due’k  <k:  Kheerraini  (1991).  The’  invisedel  analysis  is  alse) 
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applicable  to  other  vortex  profiles,  inclnding  that  of  Long  (1961),  as  stndied  by  Foster 
Dnck  (1982)  and  Foster  &  Smith  (1989). 

Little  attention  has  been  paid  to  the  stability  of  compressible  vortex  flows  of  this 
type.  On  the  other  hand  the  area  of  compressibh'  jet  flow  has  been  investigatc'd  for  some 
years  now,  the  work  of  Michalke  (1971,  1984)  being  relevant  here,  although  restricted 
to  non  swirling  flows.  Compressible  swirling  jet  flows  hav('  also  n'ct'ived  some  attf'ntion. 
the  work  of  Coleman  (1989)  should  be  mentioned,  who  stndied  the  superposition  of  a 
Rankine  vortex  on  a  top -hat  jet  velocity  fiehl.  More  recently,  Khorrami  (1991)  stndied 
a  comjjressible  swirling  axisymiiK'tric  jet,  by  assuming  tin'  incomprc’ssibh'  flow  of  Gcirtler 
(1954)  and  Loitsyanskii  (1953)  was  applicable  in  the  compressible  regime. 

In  this  paper,  we  take  cylindrical  polar  coordinate's,  {h\0,l.r).  with  the  jr  axis  lying 
along  the  axis  of  the  vortex  (which  is  taken  to  be  axisymmetric),  and  /  is  some  stream- 
wise  scale.  We  also  take  the  flow  far  from  the  vortex  centre  to  be  directed  along  the 
X  direction.  The  velocity  field  is  written  as  U*u  =  C*((/,  e,  ?('),  the  fluid  density  is 
temperature  T^T,  first  and  second  coefficients  of  viscosity  //^//,  respectiv'  ly,  and 

pressure  p^U*^p.  Here  superscript  asterisk  and  subscrijit  oo  denote  dimensional  and 
freestream  variables,  respectively,  and  also  U*  is  a  velocity  scale,  defined  below  in  (2.6), 
whilst  is  the  freestream  velocity.  We  define  the  flow  Reynolds  number 


Rr  = 


pioU:i 


(1.1) 


and  W(’  have  a  flow  Mach  numbt'r  given  by 


M  = 


u: 


1 


(1.2) 


where’  7  is  the  ratio  of  s])ecific  heats  and  R*  th<’  gas  constant.  We  also  define  the  Prandtl 
number  to  be 


(T  = 


(1-3) 


where  k*  is  the  thermal  conductivity  of  the  fluid. 
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In  tlu'  following  s('rtious  wc  rousidcr  first  tin*  hjisic  How  (st'ction  2).  which  is  shown  ‘o 
!)('  a  ])rop('r  solution  of  the  al)ov('  ('([uatioiis  of  continuity,  inomcntuni  and  ('iicrgy.  In 
section  3  W('  consich'r  th('  inviscid  linear  stability  ecjuations.  In  section  4  we  ])resent  a 
nuinher  of  nuiiKTical  n'sults.  guich'fl  hy  whi<'h.  in  sc'ction  o  we  (h'veloj)  asyni])t()tic  ic'sults 
for  larg('  wavenuinlx'rs;  throughout  this  sc'ction  W(‘  em])hasis('  a  comparison  hetwt'en  our 
nuiiH'iical  and  asyinj)totic  r('sults.  In  s('ction  G  a  ik'W  class  of  mode,  which  is  found  'o 
dc'velop  at  sufficiently  larg(’  Mach  nuinher  is  <’onsidered.  In  M'ction  7  w('  jin-sent  a  remla  r 
of  conclusions  arising  from  this  work. 
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2.  The  basic  flow 


L('t  ns  considt'r  the'  solution  ooiTrsjxniding  to  Batclu'lors'  (19C4)  similarity  solution  for  a 
swirling  wak('  flow,  ('(juivalt'nt  to  a  far  downstn'am  limit  (.;•  ;3>  1 )  of  tlu'  governing  ('([nations. 
\V('  thus  suj)i)os('  that  tlu'  solution  (•omi)ris('s  of  a  uniform  (fr('('str('am)  flow  j)lus  a  relatively 
small  [x'rturhation,  i.('. 


.s 

r  —  r.  (2.1h) 

ir  =  u\  (2.1c) 

(2.1<1) 

T=l  +  f.  (2.1(') 

l>  =  l+p.  (2. If) 


The  solution  dev(’loi)s  in  much  th<'  sanu*  way  as  tlu'  incompn’ssihh’  cas('  of  Batclx'lor  ( 19G4) 
and  it  is  found  to  h'ading  ord('r  that 


1 


"  =  V 

r  --=  0 


log(.rff( , 

S.rr* 


-Qi{v)  + 


V  0  (^  -c 


-Q2(ii) 


LV 


S.ru^'^R( 


0  •->] 


(2.3) 


tr 


C  nT}( 

2(.r?/)ir; 


(1 


(2.4) 


where 


'/  = 


r: 


S.r;/*  ^  R( 


log(.r7?f  )  4- 


Ll 


2 


Sri'*  ^R( 


{■2.6) 
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and 


(h{n)  =  e-\ 


~  ('  +  (>{v)  -  0.867]  +  2r/(?;)- 

Co  and  L  are  constants  and  Bcoc  fre<'streain  Reynolds  nninlxT  defined  by 

It(  -  * 

I'lo 


Wv  also  find  that 


r^r*  ^  rM  — 

+  2„, 2, „  . 


Sxv*JUfRf  ^  L 

whilst  the  temperature  pertnrl)ation  is  giv('n  by 


r=-  >  1  2"<4''«  >  +  S.r<'/S;r  ) 


(2.7b) 


r  ^ 


^  7^ 


8(777  2 


2^^^  (72/7  2 


+  — — c^r./-/(^/2)-8 - - - A-^c^(rf{i/2-^)),li 

i  e  -r 

+  j  J  +2- (^c^ri{(^/4) 


+  ].  (2.10) 


where  .4  and  B  are  constants,  ^  =  27/2,  ri(:)  and  crf{z)  are  tin*  expoiu'iitial  int(\c;ral  :ind 
the  ('rror  function,  respectively,  defined  in  the  normal  maniuT,  :\nd  M-y^  is  the  fna'stream 


MfU’h  number  defined  bv 


AU.  = 


u* 


u; 

Not<'  that  we  expect  A/^c  ^  A/.  Wc  also  find  it  n.seful  for  the  remainder  of  the  papt'r  to 
fh'fine  tin'  lengthscale 


i 

7-  =  7/  2  =  —  , 

^  .9 


(2.12) 


where  7',,  =  {Ax/ R/^)^  is  the  characteristic  radial  lengthscale. 

In  th('  following  section  we  consider  the  inviscid  stability  ('qnations. 


i) 


3.  The  stability  equations 

\V(’  tnkc  th('  gnirral  basic  to  be'  u  =  U(t'),  >'  —  0.  tr  =  \\  {f).T  =  To{r).j)  = 

l>[i(r)  and  consider  small  amplitnde'  pe'iturbations  to  this  How:  we-  write 

ti  =  U{r)  +  ^F{r)E  +  ()(^^  ). 

r  =  idG(r  )E  +  ). 

tr  =  \V{r)  +  dH{r)E  +  ()(d)^). 

/' =  Puir)  +  P{r)E  +  ()((^"). 

T  r„(r) +  <^r(r)£  +  0(7.''). 

P  —  pi)(i')  E  ^T(r)E  ()( e^'* ), 


wli'.'re.  •'C  1.  and 


E  ~  ('xp[/(o.r/r,>,  +  i)0  —  oe//r,,)]. 


Till'  ,ii;f)V('rning  stability  eepiations.  neglecting  the  eife'cts  of  viscosity  are'  tlii'ii 

r g'  +  /)()  ^<\  E  +  + - 1-  - ^  ~ 

p{)~pt  E  puL  Cr  =  —aP. 

2\VH  nr" 


P\)pG  +  /»()  - 


+ 


H  E  po{\v'  - 

pdf’"  +  piiG'7||(r)  —  P  +  Gr/iJ, )  =  0, 

■}EP  --  TTo  +  7p„. 


P', 

iiP 


pn 


where  we  have  written 

nir 

p  -  <t(l  --  e' )  +  . 

r 

f) 


( r  ).  f>  — 

(3.1a) 
(3.1b) 
(3.1c) 
(3.1(1) 
(3.1e) 
(3. If) 

(3.1g) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.G) 

(3.7) 

(3. Si 


and  primes  (Umote  ditfcTeiitiation  with  n'spect  to  r.  Note'  that  (3.1),  (3.2)  (3.7)  all  im¬ 

plicitly  assume  that  the  axial  scale  for  the  pertnrhtition  qnantitii's  is  considc-rahly  shortiT 
than  any  developiiK'iital  hnigthscah'  for  the  basic  flow,  as  is  the  case  if  /•„  1. 

E(piatious  (3.2)  (3.7)  may  Ix'  combined  to  yield  tin-  following  two.  first  order  <'(jua- 


tions 


(IG  ni{r\Vy  +  (\r~U'  /)(,  1 1  1  r(o-;-’  +  n‘ ) 

]  •)  CT  -f-  _  V  . 

(Ir  1-^  ')//„  /■]  />„  }-^  i(,  J 


'■  r'  'V  ^  h)  '! 


/  7  1  \  Pi) 


2t/ir  ir-^d/- 


(3.10) 


Th('S('  ('(juations  arc'  somc'what  similar  to  tho.sc'  considc'ic'd  by  Michalkc’  ( ^  .  .n  the' 

context  of  jet  Hows,  if  tlu'  swirl  vc'locity  is  nej^lc'ctc'd. 

Let  us  now  considc'r  the'  speciHc  basic  flow  of  the  trailing  line'  vortc'x.  as  discussed  in 
tlu'  prc'vious  sc'ction.  Ecpiations  (2.1).  (2.3)  (2.5)  may  bc'  snbstitntc'd  into  (3.9).  (3.10). 

and  then  if  we  assume'  |<^|  “C  jy  log.?'!,  and  using  the'  fact  that  a  sinij)!!'  transformation  <md 
invc'i'sion  of  thc'  axial  velocity  only  affc'cts  the  frecpic'iicy  of  the'  st<d)ility  analysis,  a.nd  dot's 
not  change'  the'  amplification  factor  c,  (whe'ie'  c  =  c^  +  ic,)  <md  noting  that  p.  p  anel  T  are' 
at'  orelt'r  of  O(log.r)  smaller  than  7/,  anel  ir.  then  (3.9).  (3.10)  re'duct'  to 


//(rU')'  +  or^t"  ll,. 

- - + 


o  '^r  '^  -f-  ti  ^ 


^  P. 


(3.11) 


(ir^r^)  1^,  l2r,\V 


(3.12) 


The  boundary  conditions  to  be'  applie'd  are  the'n 
P(0)  =  0.  n^O;  r'(0).  n=0. 


6'(0)  =  0.  ln|  /-  1;  G"(0).  |t/l  =  1. 


anel  G'(r)'  P{^')  bounde'd  as  r 


(3.13) 


and  the  basic  flow  may  be  taken  to  Ije 


U  =  e' 


(3.14) 


W  =  1(1  -  ), 

r 


(3.15) 


where  q  is  an  order  one  (swirl)  parameter,  and  we  have  effectively  scaled  velocities  with 
respect  to  U* . 

Equations  (3.11),  (3.12)  may  be  combined  to  eliminate  G.  yielding  the  following 
second  order  equation  for  P. 

„  „f2nlE  1  r  /\  (H'^r^)'' 

Prr  +  Pr  \  - +  -  ~  ^ - 5 -  ^ - ^ - 

r  r  r*ip  / 


+u(rH')']  I  ^  i  pnlE  W'^AP  1 ' 
(  1  r 


2nW  Wnini  (TE^r^)']' 

2  "i 

r  rV  / 


.i'2,.2y 


(ir"7' 


/I  _  +  7iirWy]  \  \2nW 

V  r  Po  )  PiS 


WHP 


PP  +7?^ 


-  -pAP  p 


{W'G-' 


(3.16) 


S(’tting  AI  =  0  clearly  reduces  (3.11),  (3.12)  and  (3.16)  to  the  incompressible  problem,  as 
considered  by  Lessen  et  al.  (1974),  Duck  &  Foster  (1980),  Leibovich  Stewartson  (1983) 
and  Duck  (1986).  It  is  also  possible  to  be  rather  more  precise  regarding  the  behaviour  of 
tlu'  sohition  as  r  — >  od;  this  takes  the  form 

'^±(\{l-P  Append,  (3.17) 

dr 

±o{\  -  PaP^P  =  0,  (3.18) 

dr 

when*  positive'  signs  ?ire  taken  for  -UKl  —  PAI^)^  }  >  0,  and  vice  versa,  to  ensure  bound- 
('dness  as  r  —>  oc. 
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Tli('  '‘order  oiu'"  proLis.  iii  requirt's  a  fully  numerical  sohitiou,  which  we  consich'r  in  tlu' 
following  s('ction,  prior  to  considering  various  asymptotic  limits  of  this  system  of  ecjuations. 
which  p('rmit  a  certain  amount  of  analytical  j)rogress. 
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4.  Numerical  results 


TIk'  sysf('in  was  treated  iiuiiierically  using  four  ditfereiit  tecliiiiciues.  The  first  was  based  on 
till'  inetliod  of  Duck  ik;  Fost(T  (19S0),  in  which  th<'  syst  'in  was  apjuoxiinated  by  second 
order  central  diti'en'iu'i's.  with  conditions  (3.17)  and  (3.18)  iinp()S('fl  at  a  liniti'  radial  value 
/■  =  '  r-  faken  sufficit'ntly  large  not  to  substantially  affect  tlu'  n'sult.  The  (h'tenninant  of 
:  1  '  y^ieui  was  tlu'ii  forced  to  xc'io  by  adjusting  the'  coniph'x  wavc'spi'i'd  c  (using  Xi'wton 
m  ra;  ion). 

Th('  si'cond  iiK'thod  was  basi'd  on  a  fourth  ordc'r  Rungt'  Kutta  schenit'.  conditions 
(3.17)  and  (3. IS)  wen'  apjiroxiinati'd  by  imposing  boundary  condition;  at  a  finite'  radial 
value'  r  =  whe'ie'  this  value'  was  again  t.ake'n  suffie'ie'iitly  large'  tei  neit  sulistantially 

affe'ct  the'  nuine'iical  re'sults.  The'  e'eunputatienis  we'ie'  jie'rfornu'el  lyv  sheieiting  the'  seilutiem 
teiwards  r  —  0  anel  we'ie'  neit  ne'ce'ssarily  e-einfine'd  tei  the'  re'al  r  axis.  The'  value'  eif  c  was 
aeljuste'f!  (again  using  Xh'wteui  ite'iatiem ),  to  e'ljsure'  the'  e'.orre'e't  be'havieiur  eif  tlie'  seilutiem 
as  r  — »  0. 

Tlie'  thirel  ine'theiel  use'el  was  base'el  em  the'  first,  l)ut  was  a  gleibal  finite'  elifh're'iie'e' 
ine'theiel.  Using  the'  afeire'Uie'ntieine'el  finite'  eliffe-n'iie'e'  se'he'nu'.  the'ii  by  de'fining  twe)  aelili- 
tieinal  epuintitie's  G  =  cG.  P  —  cP  at  e'ae'h  griel  point,  it  is  jieissilile'  tei  write'  the'  re'sulting 
se'he'ine'  in  the'  feniii 

(A-e-B)x  =  ().  (4.1) 

whe'ie' 

x=  [{C).  {C!.  {(=:■}.  )C}|  '.  14.2) 

This  sche'tne'  has  the'  eilivious  aelvantage*  eif  ge'ne'rating  4.V  I'ige'iivalue  s  (wlie'ie'  A  is  the'  nuin- 
be'r  e)f  griel  jieiints).  siinultane'emsly.  The'  principle'  elisaelvantage's  are'  (i)  it  is  neit  jieissible'  tei 
use'  (3.17)  eir  (3.18)  be'e-ause'  eif  the'  ueni  liue-arity  eif  c.  anel  e'onse'epie'iitly  Dirie'hh't  beiunelary 
e-eaielitions  we're'  applie'el  inste'ael  at  the'  emte'r  e'dge'  eif  the'  e-oinputational  eleum'dii;  aiiel  (ii) 
the-  se'he'ine'  re'eiuire's  rajiielly  ine'ie'asing  e'oinputatieiiial  re'seiure'e's  (both  in  steirage'  anel  time') 
as  .\  ine'i'e'ase's. 
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The  fourth  scIu'iik'  iiiiph'iiK'iitc'd  was  a  C’lu-hyslicv  spc'ctral  collocation  schciiic.  l)asc(l 
on  that  of  Kliorraini.  Ash  ik:  Malik  (19S9).  This  was  a  u;lol)al  incthoil.  which  generally 
iia\’<'  \(‘iy  accnrati'  ('ii!;cn\'aln('s.  hnt  also  yic  hied  a  laru;c  nniuhcr  of  sjmrions  ('n!;i’u\-aln(’s.  a 
feature  found  in  inaiiy  spc-ctral  scheiiK's.  Th<'  first  (finite  difierence)  scheme  was  (juick  and 
rohnst.  hut  hecanse  of  its  'iocal  "  uafnr«'.  iiuxh'  jninpinu:  was  often  exjx'iieiiced.  due  to  the 
fr('(inent  close  proximity  of  nei,<!;hl)onrm,u;  modes,  as  discnsseil  later  in  the  pajx-r.  The  second 
(  Rnnu;e  Kntta  )  scheme  was  also  (jnit*'  fast,  and  had  the  advantage  of  Ix-imt;  able  to  comi'nte 
neutral  and  luxir  iii-ntr.'il  modes  (and  even  st.ahle  modt's )  hy  contour  indenttition.  Ittit 
tiytnn  lx'c;ins('  it  in\'olv('d  loctd  se.archiim,  w;is  proix'  to  mode  iumpimt;.  Tin'  third  sclu'ine. 
ntmiely  tin'  ulohtil  finite  diftereiux'  scheme,  which  proved  to  he  vt-ry  rohnst.  prodnccxl  few. 
if  :my  sjtnrions  (i.e.  non  physictil)  modt's.  ;nid  <!;<'nerated  imniy  ('ietnivtilnes  simnlttmeonsly 
hectmse  of  its  xlohtil  mifuim  C'onse([n('nt ly  our  results  wt're  eenertilly  comjtntt'd  nsimt;  the 
first  sclu'tne  (if  only  ;i  limited  nnmher  of  modes  were  nxpiired)  or  the  third  scheme  if  it 
w;is  reijuiretl  to  comjtute  imniy  modes. 

We  now  [)res(  /if  a  h'w  /iiiiiierical  lesnlls  to  "ive  som<'  imlictifion  of  the  ('ffects  of 
\:ni;ition  of  cert;iin  of  the  imi)ort;int  piirtniieters.  H<T('.  ;nid  imh'ixl  in  till  our  ctilcnhitions 
W('  chose  1.4.  (/  =  O.S.  Further,  we  !>;ener;illy  found  tlnit  —  d  w;is  snffii’ieiit  with 

Aic  ^  0.03. 

Fii>;nres  1  <i.  h.  c  show  the'  v;ui;ition  of  ,u;rowth  ititi'  (oe,)  with  o  for  fix'  citst-  M  =  3. 
with  /i  =  — 1.  —2.  —3  resjx-ctively.  Tlx'se  residts  (and  ;ill  those  ])r('sented  in  this  ]);ij)er)  tire 
acenrtite  to  within  tlu'  ,i>;r:q)hic;il  /icciinicy  of  the  fiu;nres.  Bec;ins('  of  tlx'  x;re;ir  mnlfiplicity 
:ind  close  ])roximity  of  ;dl  these  modes,  we  show  tlx’  conq/nted  viihx-s  of  oe,  ;it  e;x'h  \:dne 
of  o.  f'lom  the  oiitset  it  should  he  st;ite<l  tluit  :dl  our  results  relate  to  ne'j;;iti\(-  v;ilues  of 
ir.  in  i!;eix'r;il.  we  helieve  that  instjihilities  :ue  ;dnx)sl  exclusivc-ly  confiix'il  to  this  reu;ion  of 
])ai:unetei'  sp.ace.  with  ix’rh;ii)s  ;i  few  ininoi'  ex<’eptions.  We  cert;\inly  expect  the  hiri/;est 
xiowth  r;ites  to  he  cotifiiied  to  the  ix’,ii.:itiv('  values  of  n  (this  is  .also  the  c;ise  wheti  the  flow 
is  itx'onq)ressihle ). 

We  see  in  fii/;ur<'s  1  the  followin'./;  xeneiid  ticnds;  (il  an  incre;ise  in  >j,i(iwth  lates  as  -  n 
inc/c  ises;  (ii  )  :m  ;q)]);irent  cut  ofl  v;ihx'  of  o  ;d)ove  which  no  unstahle  modes  exi-'t:  |iii)  ;i 


rc'iult'iicy  for  the  inaxiiiuun  growth  rat«‘  to  1)0  attained  at  an  inoroasing  vahu'  of  o  as  -n 
is  incr('as('d:  (iv)  many  modes  of  instahility. 

In  figures  2  a,  I),  c.  d.  ('.  growth  rates  for  tin-  eas<'  M  =  5  with  //  =  — 1.  -2.  —3.  —4.  —5 
r('s{)ecti\'ely  ar('  pr<'S('nt('d.  The  trc'iids  (i)  (iv)  descrihed  above  art'  again  ol)s('rved,  to- 
gethc'r  with  th('  rc'sidt  that  for  eorr('sj)onding  n  and  n.  the  growth  rati's  of  M  =  o  are 
substantially  r('dnc('d  rompari'd  to  4/  =  3.  Fnrtlu’r  tlu'  nppc'r  limit  of  n  of  th('  instability 
ajipears  to  be  (piite  imh'pendent  of  M.  and  section  G  confirms  this  observation.  Note-  that 
in  fignr('s  1  and  2.  the'  lower  rang)'  of  n  has  b('<'n  dc'lilx'ratc'ly  tnmcatc’d.  In  this  limit,  onr 
nmnerical  r('sidts  Ix'canu’  extrt'iiu'ly  sc'iisitive  to  grid  size.  :md  computation  rcHpiirements 
b('came  prohibitive. 

Thes('  r('snlts  sngg(’st  a  nninber  of  intt'resting  ft'atnrc's.  and  in  tin'  following  sc’ctions 
w<'  mount  a  systematic  study  for  increasing  A/,  when  ( —  )0  1. 
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5.  Asymptotic  results,  (-n)»i 
5.1  M=0(1) 

It  turns  out  that  for  this  ordor  of  Madi  nuin])er,  to  the  orders  to  which  wc  concc'rii  ourselvt's. 
the  solution  (in  particular  for  the  complex  wavespeed)  remains  unchanged  from  tlu' 
solution  of  the  incompressible  case,  as  considered  by  Lt'ibovich  ^  Stewartson  (1983)  and 
Duck  (198C).  However,  since  this  solution  forms  a  l)asis  for  the  following  subs('ctions.  w(' 
outline,  lu'iefly,  the  form  of  th('  structure  in  this  case;  full  details  can  be  found  in  the 
aforc'inentiom'd  paj)ers.  We  have  that 

o  =  ?m,  o  =  0(1).  (o.l) 

and  the  complex  wavespeed  develojis  as 

e  =  co  +  -  +  -^  +  ...  .  (5.2) 

n  |7j.|2 

Th('n  geiu'rally 

'P  —  +  9^1+  I  +  •  • .  ,  (5.3) 

|nl2 

where 

W  (  7' ) 

v5o  =  a(b^(r)-eo)+^^.  (5.4) 

r 

However  the  solution  is  found  to  be  concentrated  about  points  r  =  /’o  (critical  points) 
where  -puir  =  ro)  =  0  and  so 

o(b^(ro)-co)4--^-^  =  0.  (5.5) 

ro 

However,  it  turns  out  that  /-q  must  in  fact  also  be  a  turning  point,  and  so 

V^'(ro)  =  0.  (5.G) 

i.e. 

(^U'(rn)+(~)  =0,  (5.7) 

V  r  /  () 
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(wIk'i'c  a  subscript  zt'ro  hen'  and  li(T('aft('r  denotes  evaluation  at  r  =  Co  ).  Ecpiation  (a. 7) 
tin'll  serves  to  d('t('rmin<'  rn.  and  Ik'iici'  o,  may  he  deti'iniiiu'd  from  (5.5).  The  key  l<'n,c;th- 
scah'  inside’  the  critical  lave'r  is  a;iv('n  hv 


;?  =  (/■- 


'■(I  hd  ‘. 


and  till'll  on  this  scale' 


sigiiln 


Y  r  I 


{— ncv  +  }  +  ()( 


will'll' 


r'l  —  — '■>^'1  • 


:5.io) 


ro('n  )  -  — 


(5.11) 


For  consistency  it  was  shown  hy  Leihovich  Sti’wartson  (19S3)  and  Duck  (19SG)  that 


,  _  (ir^r-’)' 


(an;,  + 

( 1  +  )cj^n- 


(0.12) 


1  ai  fhi'i  .  on  the  7?  0(1)  scale  the'  eigenfunctions  scale  as 


P  P  +  ...  . 


1 5.13a  1 


G  uG 


(5.131)1 


and  then  P{B)  is  descrihi’d  hy 


Phh-{2\'^ 

l  Lci 


l!(n,  )B'^  ]  sign(  n  )( 1  +  (\  ~rl )  t  — 


2(A,)^-' 


(5.15) 
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Eciiiatioii  (5.15)  has  a  sohition  which  may  l)e  writton  in  terms  of  Wc])cr  parabolic  cylimh'r 
functions.  and  so  if  we  demand  the  solution  decays  as  |^|  — >  oc.  then  m  must  he 

an  integ('r.  yielding  the  following  r{'sult  for  c^: 


Figure'  3  shows  the'  variation  of  growth  rates  for  M  =  3.  for  n  =  —1  (tlu'  h'ast  unstable 
mode'  shown)  up  to  =  —10  (the  me)st  unstable  me)ele'  she)wn)  as  ce)mpute'el  fen-  the'  fidl 
syste'in:  in  all  case's  the-  most  unstable  nmde  for  e'ach  e'aliie'  of  it  is  shown.  Using  tin- 
asymi)t(5tic  residts  above,  we  she)w  the  corresponeling  re'sidts  in  figure'  4  (in  particular  we' 
se't  in  =  0  in  (5.17)).  The-  comparison  be'twe'e'ii  figure's  3  anel  4  re've'als  ge)e)d  agre'enie'iit  in 
the'  gre)wfh  rate’s  at  large'r  value's  of  — n.  altlienigh  at  smalle'r  value's  of  — n.  comparison  is 
le’ss  goenl  than  the’  comparise)n  e)f  incompre’ssible’  mnne’iical  anel  asym])te)tic  re’sults.  This 
le’ss  satisfactory  agreeme’iit  may  be  attribute’el  dire’ctly  to  the  e’ffe’cts  e)f  comi)re'ssibility.  but 
for  [t?!  sufficiently  large,  the  asymptotic  results  pre’sente’d  in  this  sub  se'ctie)n  are’  edtimate'ly 
ap])re)ae'he’el. 

In  the  following  sub  section  we  be’gin  to  inemrpeirate  cennpre’ssibility  inte)  e)ur  asymp- 
te)tic  eh’seription. 


5.2  M=0(ln|i) 

The  ke'y  enpiatieni  (5.14)  was  eibtaine’d  l)y  taking  the’  ()(|7/|^  )  terms  in  eepiatieni  (3.1G). 
If  |7/|  ^  1.  anel  we  ne)W  pe’rmit  M  to  grenv  in  magnitude’,  adelitional  terms  will  ultimate’ly 
e’lite’r  (5.14)  wlie’ii 


i\/  =  l7/|UV/,  A/  =  0(1). 


(5.1S) 
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wlu'ii  th('  coefficient  of  P  will  include  the  effects  of  coini)ressil)ility.  The  inoflific'd  ('(|nation 
is  tin'll 


Prr  -  sign(n)|2|^ 


y'( 


2d  Cl 


+ 


(rir);  +  or;ir'i 


,.10, 


All  otlu'r  (luantities  (spi'cifically  co  and  C)  )  remain  unchanged  from  thos*'  evaluati'd  jiri'vi- 
ously.  This  equation  may  be  transformed  to  the  same  form  as  (5.15).  viz 


+  ( 


AqC^  +  A2) 
2(A,)^ 


(5.20) 


wlu'n'  A()  and  A]  are  givc-n  hy  (5.1G),  whilst 


A2  =  .sign(  n ; 


(rir),',  +«r5c,; 


orjr, 


'■() 


Consequently,  using  our  previous  arguments 

( Ai )  2  (1  +  2m )  +  A2 


('2 


Af 


m  =  0.1.2, 


(5.21) 


(5.22)) 


Figure’  5  shows  a  comjiarison  of  “exact"  growth  rates  (obtained  using  the  numi'rical  ap¬ 
proach  of  si’ction  4),  shown  as  a  broken  line,  with  the  asymptotic  results  of  this  subsi'ction: 
h('re.  the’  cheise’ii  value’s  of  n  are  -5,  -10  and  -15  anel  we  set  M  =  |n|7  i.(’.  A/  =  1  in  e’ach 
e-ase-.  The’  agre’enient  is  sati.sfactory,  and  inijiroves  as  n  increase’s.  Figure  6  sheiws  the’ 
e-orre’spemeling  comparison  for  rv;  in  this  case  the  agreement  is  e’xcellent. 


5.3  M=0(|n|t) 

As  A/  (anel  hene-e’  A/)  increase  in  magnitude,  the’  “A2  e'ompone’nt"  of  ('2  increase’s  (as 
\P  )  and  be’cenue’s  sufficie’iitly  larger  than  “Ai  cennpone’iit” .  Simultaneemsly,  the’  ceietfiede'nt 
eif  P,  in  (3.1G)  will  grow,  anel  furthe’r  te’rms  in  the  coe’fficie’iit  eif  P  will  be'e'emie  significant. 
The’  ne’xt  impeirtant  re’gime  for  M  is  whe’ii 

A/  =  |n|^A7,  XI  =  (){\).  (5.23) 


IG 


Although  tlu'  k(’y  r.f.lial  srah'  n'inains  R  =  0(1)  (see  (5.S)).  the  sc'ric's  dovolopnu’iit  of  r 
and  ^{R)  is  now  alt<'r<'d,  and  is  instead 


I  '  '  I 

(  -  f'o  H - f  7^ 

7/U 


(5.24) 


and  therefore. 


,  y2  ,  , 

^  +  ••• 

77  1  77  2 


-fV7l 


n  -1 


-n7''2sign(77) 


^|-oe3sign(77)  +  77o(n))-R'^}  + 
77  2  *'  Z  ) 


(5.25) 


Considc'ring  {)(|77l-'  )  ([uantities  in  (3.1G)  requirt's  that  the  sum  of  these'  t('nns  is  ze'ro.  and 
so  this  h'ads  to. 


7-2  =  -sign(77)— — 
2o7-o 


1  +  d'^'o 


(5.20 


The  ('(piation  for  P(R)  is  again  ol)tainefl  by  taking  the  0(|77|2  )  tc'nns  in  (3.16).  n.'unely 


fn 


(5.27 


Using  the  standard  transformation 


P(i?)  =  P  (P)exp 


1  APW^iR 


then 


^*Rli  -  {2  sign(77)^^ 


‘do(d)  \  /  1  + 


2oe 


)C-^) 


(5.29) 
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or 


whcrc 


A:, 


1  ^2  I  A()r2  +  A.t 

2(A,)^ 


Al^ 

2/-„ 


P  =0. 


(5.30) 


:5.31 


Since  (5.30)  is  ;iu;;un  ;i  form  of  W'elx'r's  eqmitioii.  then 


C', 


( A|  )  -  (  1  -f  2>n  )  +  A.t 


}■ 


,»  =  0. 1.-2. 


15.32) 


Fii2:nr('  7  shows  ;i  comi):irison  of  nsynipfotic  results  (solid  lini's)  with  nninerical  i (-suits 
(hroken  liiu-).  for  n  =-5.  -10.  -15  with  M  =  ln|~  (i.e.  M  =  1).  .Ayain  tin-  coin])arison 
l)ecoines  markedly  better  as  |/(|  increases.  The  corn-spondinti;  distiibntion  of  e,  is  a^aiii 
Sood.  and  shown  in  fi£>;ure  S. 

5.4  M:=0(|n|^) 

In  this  case  it  is  (jnite  clear  from  tin-  jnevious  subsc-ction  tlnit  when  M  =  ()(|//|s  ).  then  tin- 
cj  term  (see  (5.2G)).  which  .throws  as  will  Ix'come  comparable  in  mas:nitnde  to  tin-  t('rm 
involviii'j;  C)  (which  is  imh-jx-ndent  of  M ):  in  addition  to  this,  since-  f  {  throws  as  A/'  as  AI 
increase's,  then  this  term  will  also  ix-come  comparable-  to  the-  C|  te-rin  whe-n  M  =  ()(|n|7  ). 
In  this  case-,  it  is  straip;htforward  to  show  that  /?  =  ()(!)  re-mains  the-  aj)pre)j)riate'  raelial 
scale-,  whilst  the-  wavespee-el  e'X])ands  as  fe)lle)ws 


implyin<>; 


e-,  ('2 

•(.  +  -  +  y-TT  +  .  .  .  . 

bd^ 


.  I  T"  , 

r  I  +  rrr  + 


(5.33) 


(5.34) 


P)  2 


If  We  wiite 


M 


KM.  M  =  (){]}. 


(5.35) 


IS 


rlu'ii  the  ('(luatiou  for  P{R)  is 


P 


HR  ~  i  I" 


/•If  )||  +  o /'i'; f  i'  \  /  2f f  (I 


'  II  r  I 
1  +  '  I  ■  /■,■; 


—  siu;u(  II ) - 

orl  'n 


m.. 

^  'nr  I  ^ 


I'd  ' 


'  II  r  I 


+ 


f  sii!;n(  n  ] 


'  II  r  1 

/(/•If 

Cl  + 

'>'dfll  \  ( 

.u^if,r\ 

1 

'  II  r 

> 

/•i.  / 

+2| 

/  1  +  o  “  /■,•; 

^  'nr'l 

'lir- 

r  I 


=  0. 


Abilin  iisiiie;  a  standard  f raiisfonnation  of  thr  form 

P{P]  =  P(P)cx\) 


1  \n\-<M-\V-R 


I'u 


li'ads  to  till'  ditfcrontial  <'<iuation 


p 


HR 


f  iV  /'o^'l  ^'^'nrl  '  '-M  ’ 


+ 


'  -',^-1  A-'  )  +  V  2n,  ) 


_  9  ^  ^0  d-  ''pf  n  ‘  y  2 


/’d  r  1 


)■ 


r  I 


+siRn(  // )  ( fl 

'  '"d  r  I  ^  ^  '’<1  ^  r  I 


/  1  +  '*''n  ^ 

^  'dr.  ^ 

^  '-dd. 

P  =  0. 


( 5.30 ) 


( 0.37 


(5.3S) 


If  our  assumption  n'l^ardintz;  the  importanci'  of  thr  /?  =  0(1)  M-alc  is  to  Ix'  i-onsistcnt . 
( altliotif>;li  SIX’  the  comnH'nts  Ixdow  roG;ardin,i;,  R  —  ()(|n|^))  tlirn  the  corflii-irnt  involving; 
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and  usinu;  j)rcvi()ns  ai7!;uniciits.  W('  must  liavo 


(  A|  )  2  ( 1  +  ) 

A() 


m  =  0.1,2... 


(.3.42) 


.\otc  that  alth()Ut>;h  tlu'  transfonuatioii  (5.37)  suggests  growth  as  7?  — >  oc.  this  is  inon'  than 
offset  i)y  the  decay  of  tlie  jjaraholic  cylinder  functions,  allx'it  on  a  hargt'r  h'ligthscah'  7?  = 
()(|nj^  ).  Xot('  too  that  tlu'  transformation  (5.37)  is  <'onsist('nt  with  that  ust'd  ])r('viously. 
namely  (5.28). 

Figur('  0  shows  comparison  he’tween  the  above  asymptotic  rt'sults  (solid  lin<')  and  tin* 
"exact"  iiunH'rical  results  (broken  line),  for  tin*  most  unstable  growth  rate's  wln'ii  n  —  -5. 
10.  15  (.17  =  |n|K  i.e.  .17  =  1  in  all  cases).  Again  the'  agre'enu'iit  is  se'e'ii  to  improve  as 
-11  increases.  Figure'  10  shenvs  a  e-ennparison  be'twe'e'ii  the  ce)rre's])e)neling  e  ,  s.  fe)r  the'  abewe' 
case's  anel  inelie'ate's  ge)e)el  agre'e'ine'iit  be'twe'e'ii  emr  asymjiteitic  and  nume'rie-al  re'sults. 
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It  turns  out,  liowrvc'r.  that  this  oifh'r  of  M  marks  a  \vat('rsh(‘(l.  WIk'ii  M  is  not 
large,  lioth  roots  of  r|  must  Ix'  comph'x  (an<l  from  our  previous  discussions,  the  How  is 
unstahl(').  As  M  increast's.  how('V('r,  (wentually  t«Tnis  insid('  tlu'  squan'  root  term  in  (3.39) 
will  ('Vt'iitually  Ix-conu'  positive',  and  h('nc<'  the  two  roots  of  C|  will  cease  to  he  <-omi)lex 
conjugate  pairs,  hut  will  hoth  Ix'conu'  n'al.  In  particular,  this  will  occur  wIk'u 


(3.43) 


To  illustrate'  the  stahilisation  of  the'se'  mode's,  in  Hgure's  11  a.  h.  e-.  we-  she)w  the’  variatiein 
e)f  the’  gre)wth  rate’s  with  M  for  the’  e-ase-s  ii  =  —3  anel  o  =  2.5.  ii  =  —10  ami  e>  =  3. 
>1  =  —13  anel  ee  =  7.5  re'sj)e'ctive’ly.  We  ede'arly  se'e’  the'se’  me)ele's  he’e‘e)niing  iie’utrally  stable’ 
at  Hnite’  value’s  e)f  M.  Fe)r  ceniiparise)!!,  in  figure’s  12a  12c  (e’e)rre'sj)e)nehng  te)  Hgure’s  11a 

lie-  re'sj)e'e’tive'ly)  we'  she)w  the’  cen're'speuieling  re'sults  freun  se)lutie)ns  e)f  the'  full  syste’in  (3.11 ). 
(3.12).  At  the'  lowen'  value's  e)f  d/.  the're'  is  ge)e)el  e-euTe'latiem  he’twe'e'ii  the'  twe)  se’ts  e)f  re’sults. 
He)we’ve’r.  as  the’  Mae'h  numhe'r  incre'ase's.  anel  as  the'  eu'ele'r  e)f  the'  me)ele's  (i.e’.  tii)  im-rease’s. 
the'  cen'ie'latie)!!  ele’te'rie)rate's.  W'e’  cenisieh'r  this  latte'r  j)e)int  first  hy  e’xamining  the’  he'havie)ur 
e)f  the’  highe’r  e)rele'r  mexle’s. 


5.5  The  structure  when  m  =0(|n|?)  ,  (M=0( |n| ^  )) 

One’  eh’tail  that  has  ne)t  he'e'ii  ce)nsiele’re’el  se)  far.  is  the’  he’havie)ur  e)f  the’  meeele's  as  iii  ine’ie'ase's. 
This  feature  was  cemsidere'd  hy  Duck  (198G).  anel  we  cemsieh'r  this  as])e’ct  he’re’.  fe)r  the’ 
impe)rtant  case  wlie’H  M  —  0(|/e|2).  c  ele've'le)ps  in  the  same  manne’r  as  in  the'  pre'vie)us 
suh  se'e'tie)!!.  anel  ele)se  te)  the  critical  laye'r  y  take's  the’  form 

1 

y  —  yl  +  ■ — +  •  •  ■  •  (3.44) 

|n  1 2 

e'o  re’inains  unchange'el  fre)m  that  give’ii  pre'vieuisly.  he)we've'r  eq  will  eliffe’r.  In  this  e-ase’,  we' 
e-e)nsiele'r  the'  le’ngthscale 


/?  =  ('■-  =  0(1). 


(3.43) 
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and  so  oil  this  scale 


r'l  —  —^'<'1  H - rol'i*'-'*  • 


=  —  0C2  + 


- . „)R\ 


3! 


r  0 


Snlistitutioii  of  these  ('xpaiisioiis  into  oiir  ,!>;oveniiiig  eipiatioii  ,e;iv<’s.  to  lea(liii«' 
eii;;eiifutictioii  ('([nation  of  th('  forin 


Usiinr  the  standard  transfonnation 


P(R)  =  P*[B)f.rp 


HI 


q\{R)(lR 


wv  obtain 


Ho\V('V('r.  sinc<’  |/)1  is  larg(',  (5.50)  may  be  ai)i)roxiinated  liy 


=  0. 


\vh('re 


'  ()y  1 


4roH  o^  +  ( 1  4-  )(b)'r’  i  “ 


-sign(n)2Tr,fr,^d/^ypi  + 


-i  ^2  ^ 
0  r  I '  0 


Thus  (5.50)  may  he  written  as 


Pi  + 


^  ( 1  +  ^'^  )  ( 1  +  ^'^  ) 


P*  =  0, 


(5.4G) 

( 5.47 ) 

ord('r,  an 

(5.4S) 

( 5.49 ) 

(5.50) 

(5.51 ) 

(5.52) 

(5.53) 


90 


with 


2[4roir|f  -(1  +  (v^-o)(n’'^?--)[)]sigu(//) 
-2[(1  +  +  A/'U',)V4]nrisig;n(?0 

'■ov^o  ('■<)) 

r-v"(n))sigii(//)l  "  ~ 


( 5.54 1 


(5.55) 


( 5.5C ) 


(5.57) 


P  =  \n 


(5.5S) 


Etiuation  (5.53)  i.s  now  in  a  form  suital)U'  for  a  WKB.T  tyjx'  of  approximation.  This  ('cpiation 
has  four  turning  points,  howevx'r  for  the  rango  of  ^  n-fpiinxl  that  W('  lu'c'd  consich-r  only 
tho  turning  j)oint,s  at 


2u  2)1  P  ^  I' 


1 


1 

2 


(5.59) 


For  large  h  the  WKB.I  solution  is  given  l)y 


=  q 


xp/(, 


.4i  exp  /  {thpdp  4-  .4-2  exp  /  —{ilipcjp 


/ 


(5.60) 


where 


l< _ 


(5.61) 


The  treatnu'nt  near  the  turning  points  is  standard  (see  for  example  Duck  1986)  and  leads 
to  the  following  dispersion  relationship  for  rj 


//,  =  -~(1  +  27n),  77/  =  0.1.2,... 

4'  ' 


(5.62) 
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where 


1 

~-4‘^rT^7yi7T' 


1 

2 


(5.G4) 

( 5.05 ) 

(5.GG) 

(5.G7) 


(5.GS) 


Here.  F(~.k).  E(~.k).  n(^.o^,A')  denot*'  roniph'te  elli|)tic  iiitf'gral.s  of  the  fir.'^t.  seeoiid 
and  third  kinds  resjx'ctively. 

The  system  (5.G2)  was  solved  using  Newton  iteration,  and  n'sults  for  variations  of 
growth  rate  witli  M  w('re  computed  for  tn  =  0.1, 2, 3, 4. 5  for  the  cast's  ti  =  — 5,o  =  2.5 
(tigure  13a),  n  -  10,  o  =  5  (figure  131)),  ii  =  -15.0  =  7.5  (figure  13c).  The  modifii'd 
theory  of  this  snl)section  does  indicate  some  important  imp.rovemeut  in  tin'  coni|)arison 
with  figures  12;  in  ])articular  the  "Inmching  iip”  of  tlu'  moch's  with  an  incn'asc'  in  ordt-r  is 
captuK'd.  However,  although  for  n  —15  there  is  <j\iit<'  good  coriH'latiou  for  the  lowt'r  ordi'r 
modes,  as  ni  incrc'ases,  the  various  wiggles  ohs('rved  in  hgun's  12  are  not  described,  and 
more  im|)ortantly  onr  asymptotic  res-ilts  do  not  capture  tlu'  instaoility  shown  in  figures  12 
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Ix'VoikI  al)()ut  M  =  S.  Iiulcrd,  figun*  12  indicates  instability  as  M  iiicn'asc's.  \V('  thcn-foic 
coiicludc  that  tlu'  iiatnrt'  of  those'  mode's,  as  M  —y  oc  is  se)nie'\vhat  elitfe're'iit .  It  turns  e)ut 
fre)ni  the'  analysis  eef  the'  fe)lle)\vin,G;  se'ctie)n  that  as  d/  — >  dc.  a  furtheT  class  eef  meeele'  e'lncrse's. 
ejuite'  elistine't  freun  the)se'  e-e)nsiele're'el  se)  far. 


6.  Centre  Modes,  M=0(|n|) 

Hrn’  \V('  rxainiiu'  iiiocU's  i'xhi])it('(l  whc'ii  the  Mach  iiuinhcr  is  of  ordtT  n.  i.('. 

3/  =  |;;|Ir.  If  =  0(1).  (C.l) 


As  not('(l  pn'viously.  as  M  oo,  the  fcattiros  of  the  iiuuh's  exhibited  at  lower  Mach 
nuinlxTs  C('as('  to  I'xist.  In  fact  the  analysis  of  snh  sc'ction  5.5  sne;g('sts  that  tlu'  inodi's 
in  this  n'ttiiiK'  nndt'r  {)res('nt  consid<'iation  an'  neutral.  Howev('r.  our  nuiiK'iics  bi'lie  this 
and  point  towards  th('  ('xist<'nc('  of  c.c.ntrr  viodc.-i.  by  which  wt'  nu'an.  th('  eigi'iivahu's  are 
d('t('rniin('d  primarily  by  conditions  close'  to  the  axis  of  tlu'  vortex,  r  =  0.  in  a  nnuiiu'r 
simih'ir  to  that  of  Ste'wartson  and  Brown  (1984,  1985).  He'n'.  tlu'  coinph'x  wavesjx'e'd 
(h'velops  as 

r  =  ("(I  +  —  +  .  .  .  .  (6.2) 

tr 


(;ni  e.xpansion  that  can  be'  ve'iifie'el  ti  po-Hcriori)  and  the'ie'fore' 


y  I 

=  //yd  + - 1-  .  .  . 


=  /t|o(r  -  ((i)  +  ^}  +  } 


(6.3) 


In  e)rele’r  te)  be'  e-e)nsiste'nt  we'  must  luive'  that  =  6)  =  0,  i.e'. 

Co  =  1  +  -.  (6.4) 

We  must  ne)W  ge)  on  te)  finel  C|  in  enir  e'e)mj)le’X  wave'spe'e'el  e'xpansie)n.  The-  se)lutie)ns  are', 
in  the'  main.  trapi)e'el  in  the'  re'gie)n  arenmel  r  =  0,  but  we'  be'gin  by  she)wing  he)W  the'se- 
se)lutie)ns  are'  e'e)nne'e-te'el  te)  the'  e)ute'r  He)W.  The'  fle)w  is  eliviele'el  inte)  the'  fe)ur  re'gie)ns.  whie'h 
are  e'e)nsiele'reel  in  turn. 

We  first  be'gin  e)ur  analysis  at  the'  e)ute'r  n'gie)n  e)f  the'  How.  wlie'ie'  c  =  fl(l)-  ^He' 
ge)verning  e'epiatie)!!  h;is  the'  fe)rm 

Prr  ~  +  \n\^Jl\^lP  =  0. 

/’ 


(6.5) 


Using  a  transforiiiation  of  the  form 


P  =  P*(r)cxp-\  I  - ; - dr 


’  ro 


ir 


and  writing 


\()  —  —^P^l^  + 


WPM 

4r2 


r('dnc('s  (G.5)  to 


p:,.  -  \v\\oP*  =  0. 


(G.G) 


(6.7) 


(6.S) 


We  dc'fine  the  critical  point  t'a  l)y  \n(7‘o)  =  0  and  wlien  r  >  ;'o.  \o  <  0  and  the  approximate 
solution  to  (6.S)  is  given  i)y 


P*  = 


(6.9) 


where  Ei  is  an  unknown  constant. 

We  have'  also  assumed  that  [— \o]5  is  slightly  imaginary,  with  its  imaginary  part 
positive'.  This  condition  is  to  ('iisun'  that  th<'  solution  is  houneh'd  as  r  — »  dc.  oth('rwis('  if 
/7n{[  — \()]?  }  <  0  W(’  n’place  by  '  —  in  (6.9).  Looking  at  this  i)oint  from  a  more'  physical 
perspective,  we  require'  waves  to  propagate  out  from  the'  critical  laye'r  and  not  in  toward 
it. 

In  orde'r  to  lie  able  to  match  (6.9)  with  the  solution  in  the  ;•  <  f’o  rt'gion.  wc'  must  first 
('xainiiK'  th('  solution  of  the'  fleiw  in  the'  transitiein  laye'r.  areiunel  ?■  —  r'o.  The  k'ngthscale  in 
this  re'giein  turns  out  tei  be 


r  -  I'n  =  i?|/7|-U  i?  =  0(l), 


(6.10) 


anel  use'  eif  this  scaling  re'eluces  emr  geive'ruing  e-eiuatiein  tei  the'  Airy  ('epiatiein,  which  has 
seilutieins  eif  the'  feirin 


P*  =  Fi  Ai(\,',(7-o))^.ff)  +  F2Bi(\,',(7'n))^'^). 


(6.11) 


Standard  analysis  and  inatchine;  with  (G.9)  yields  the  n-sult 


F,  =  ,F,. 


(G.12) 


In  th('  rc'gion  r  <  ro,  Xo  >0  and  lu'nec'  the  WKBJ  solution  lu'if'  is  givtni  hy 


■xp{-  f  \„\‘ 

J  rn 


I" I  - voj 


(G.13) 


As  w(>  tipj)roach  th('  transition  la\'('r  this  solution  iimst  niatrh  with  (G.ll)  iinplyino' 


D]  ==  —iD^  \/3. 


(G.14) 


Note  that  as  r  — >  0 


/F\/^•■■^hd'^ 


F\/'rFtd- 


J  ^  ^  -u  ^  ^  f  qFM  r 

i  <  D\Z]  ex^)  y - - j  +  D2 Z2  <'xp  - j 


(G.15) 


wh('re. 


(G.IG) 


We  now  eonsid('r  tlu’  h'ligthscale.  r  =  ()(|?/|“i  ).  by  s('tting 


r  =  r|n|^  =  ()(1), 


(G.IS) 


The  fxoverning  equation  in  this  region,  to  letiding  ordt'r.  is  givt-n  by 


Pr,  ~  \r>\q'M  rP,  + 


42t/FU  1 


^ - -  F  =  0. 

rt  r- ^ 


(G.19) 
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where 


V'’i  =  :^y’o(0)^■^ 


(6.20) 


and  as  previously,  it  is  necessary  to  use  a  transformation,  which  is  of  the  standard  type. 


P  =  P**exp{^  I \v\,frMd7-Y 


(6.21) 


which  h’ads  to  the  final  form  of  the  governing  equation  in  this  n'gion.  viz 


P**  -p  - l_  _  q  r  M  ^  ^ 


(6.22) 


This  ecpiations  yields  solutions  of  the  form 


(6.23) 


wlu're  /  and  K  are  modified  Bessel  functions,  written  in  standard  form,  with 


2  V 


(6.24) 


Now,  as  r  —7  oo 


C,;- 


q^M  |n|7-‘^ 


7.-2  ,  — r 

2nq^M  |r).|r'^/4 


q^M  1/(|/- 


+C2rY  —=^ - 


).  (6.2i 


and  comparing  this  with  ('quation  (6.15),  we  see  that  in  order  for  these'  solutions  to  match 
and  taking  into  account  equation  (6.14),  then 


C,  =  -ins/3Z'fC2. 


(6.26) 
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Also.  \v('  tak('  iK)t('  of  the  fact  that  as  r  — >  0.  then 


■2\// 


Ci 

r(n  +  i)  '  ^ 


yy 


(G.27 


\vh('  n'  Ji  =  M~\}i\/4  and  r(;-')  is  the  gainina  fnnetion.  <lefined  in  the  standard  fashion. 

X('xt  we  set 


H  =  =  0(1). 


(G.2S) 


and  in  this  regime  ^  (hwidops  as 


y  — 


y  I 
// 


+ 


(G.29) 


with  th('  governing  ('qnation 


Plif}  +  |r;l‘ 


wliich  has  th('  solution 


4(\^q‘‘  +  dory 


+ 


2q'\\r 


-oe,  +  iy''(0)/?2]-^  ■  [-oe,  +  ^y;;(0)7?-’] 


1 

7?-'J 


P  =  0. 


P  =  - ^^-7— F(-,;-2:|"I  + 

(1  -  Tn^y‘  '  '  "  '  ' 


(G.3()) 


(G.31) 


wh<T('.  F  is  tlu'  hyp(’rgeoni('tric  function  in  standard  notation,  and 
^  _  yyo(O) 

OC| 


^  1 


:^U"I  - 


•^2  —  k'l  +  -s)  —  //, 


•S  —  (  |//  1'^  +  . 

r/T 

/-/(/-/  A  1)  . 


(G.32a) 

(G.321)) 

(G.32c) 

(G.32d) 

(G.32e) 
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with 


i  4(yo  +  4(j~(\ 1 

I  \{r[:m7' r 


I  G.33:i ) 


(G.331)) 


Xoti’  til, -it  this  ('(Illation  is  v('iy  similar  in  form  to  that  found  hy  Stcwartson  Caja'll 
(1985).  This  solution  is  of  (•ours('  finit('  at  R  ~  liy  tlu'  (h-finition  of  th('  hv])('ri;('omcti ic 
function,  and  therefore  satisfies  our  lioundary  conditions  at  the  centre  of  the  \'orfex.  As 
7?  liecomes  larG;e 


P  ~  R^A^R^  +  AiR-^). 


(G.34) 


imitcliinu;  this  witli  (6.27)  and  makin.s;  use  of  tC.2G)  we  se('  that. 


r(.s/2  +  l)n.s/2)|n 
(7.^1 


.-ii 


-/;r.4.2v^. 


(G.35) 


whicli  u;ives  us  a  r('l;itionship  l)etw(>eu  .4i  and  .-li-  We  now  pjo  on  to  d('t('rmine  ci  explicitly. 

We  hegin.  hy  writing  P.  as  R  cc.  in  terms  of  its  asymiitotic  (-xpansion  in  flii" 
limit,  i.e. 


P 


(T/?-’)^  +  ^ 
(1  -  TRA" 


S(l( 


-T/?")"-'^  + 


(G.3G) 


where. 


(G.37a) 


(  G .  i?  ( 1 )  I 
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and  and  an-  both  power  series  expansions  in  1/T/?',  with  their  h'ading  order  tennis 
unify.  Matehini;;  tlu'  two  huge’  R  solutions.  (C.34)  aiul  (G.3C)  yields  the  relationship. 


yip 

.d. 


Sd 


T)\ 


(G.3S) 


Ih)W('ver.  since  l-'^l  0  then  |T''|  0.  which  in  turn  iinjilies  tlnat  |>.|)|  j(\i|.  and  tlu'refore 

to  leadiini;  order,  sn  =  0.  Hence.  w<'  see*  that.  :7(|n[  +  .s  j  —  fi  must  Ix'  ('itlu'r  a  negative' 
inte'ge'i'  eir  ze'ro.  i.e. 


fi\>  =  +  -V. 


(A’  >  0). 


(C.39) 


for  infe'ger  .V.  wlie're  /hi  elf'iiote's  a  first  orde-r  approximation  to  /;. 

Sine'e'  .s  is  imaginary,  the'  ahewe*  e'xpre'ssion  suhstitute'el  into  (G.32i')  and  (G.32a)  give's 
the'  le'aeliii'g  orde'r  imaginary  e-ontriliution  to  C].  i.e'.  it  is  ohtaine'el  from 


oOn(O) 
Se'>//(|(  //()  +  1  ) 


(G.4()) 


It  is  alse)  peissihle'  to  obtain  the'  e'orre'e'tion  to  this  te-nn  by  re'taining  the'  ne'Xt  eirele'r  a])pro.x- 
innitiein  te)  ji.  On  se'ttin<g.  fi  =  /hi  +  A//,  when'  1A//|  'C  1.  the'ii  we'  obtain 

(  -  l)'^'+‘|7h!.V!(.s  -  1)! 


(s  -f  |nl  +  .V)! 


A/h 


(G.41a] 


^1 


(-l)''|n|!(.s  +  .V)! 


(G.41b) 


r'lnel  tlie'ii  (G.3S)  givc's 


''  .4,  .s!(,s-l)!.V!(in|4-A’)! 


(0.42) 


Sine'i' 


//(//  f  1)  = 


4etC| 


(G.43) 
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and  thcifforc.  if  \vv  allow  //  =  //n  +  A//  in  tin'  above,  we  see  that 


dO) 


'■|  = 


1) 


S(»//(i( /'ll  f  1  ) 


d/'ii(  /'ll  h  1  ) 


A/' 


((i.44i 


Therefore,  we  m;iy  now  write  e|  ;is 


'1 


^n(O) 


So//|,(/'/|)  4-  1 ) 


1  4 


.4, 


( 


Aj  \  4/'(i(//i)  4  1 ) 


I  >  a  |"|  +  A'  )!  !  >  -  .V  l! 
s'l.r^Ti!  A!oe  A  Yii 


4/',i(  /'ll  1  I J 

where,  till  varitthles  in  this  ('xpression  htive  tx'c'ii  previonsly  ilefined,  Howevi'i .  the  eorn'i  tion 
term  IS  found  to  lx-  of  severtd  ordc'rs  of  m;i,ii;nitnde  snitdler  th:m  the  leadimi  onlei  t.-rm  of 
'1  ■ 

In  figun' 14.1  we  show  th<' varititioii  of  the  ,<;rowth  rate  with  o  foi  t  he  rtise  n  --  -0.  M  - 
5.  This  is  fo  he  eoniptired  with  figaire  2e.  eompnted  from  the  foil  sy'^tem:  l•oml)<■lrison 

hetwei'ii  the  iisymptotie  ;md  fully  nnmeri«':d  rc'stdts  is  snr])risin,i;ly  liood.  considei imr  the 
simtlliK'ss  of  — n  ;aid  M.  Fifftin's  141)  shows  the  ivrowth  r.ate  vtniation.  jnedieted  from  the 
analysis  of  this  snhseetioji,  for  the  <")se  ft  =  — ]().  M  =  10  and  is  to  Ix'  cotiiptiied  with  tlje 
fully  immerical  restdts  shown  in  H/i;nre  15.  .•X.t^ain.  the  results  art'  ('neonr;i,i;in'j:,  in  ]);ii  timhn 
till'  m;i,i>;nitnd('  and  loiaation  of  the  imiximnm  .urowfh  rtite  is  (piiti'  iircurtitely  ptedieted. 
h  nfortnmitely.  onr  att('mj)ts  ,at  ;i  nnmerietil  solution  of  (.1.11).  (3.12)  with  ;i  reasomihle 
nmnlx'r  of  i;rid  points  for  tlu'  cttsi'  n  =  —15.  M  -  :  15  ])roved  unsueeessfid.  due  to  a  harye 
numher  of  spurious  (i.e.  highly  c,rid  (h'pendent )  modes.  Indeed,  it  was  found  thtit  ;is  M 
w;is  incretised.  our  numerical  srheiiK’  often  ftiiled  jdutiptly.  with  ;i  sudden  generation  of 
hire;''  numhers  of  s])urious  modes. 

Referrinc,  a,e;;iin  to  fi,<;tir('s  14.  we  see  tlnit  the  ui)J)it  neutr;il  j'oint  is  clearly  '.ee]!.  This 
occurs  ;it  o  —  n/'/.  th<’  point  ;it  which  "  —  0  (see  (G.lli).  The  fidly  tiumerical  residts 
throu/2;hout  this  pajx'r  all  vi'iy  clearly  show  this  icsult.  which  apjieais  to  .ajii'ly  uni\ersally. 
as  in  th('  incompressihle  case. 


Adilif ioually  14  show  a  distinct  lowin'  ncutrtil  point.  This  corresponds  to  the 

\ahie  of  o  where  >-'”((1)  =  0  innnely  o  =  — <y/2.  At  this  location.  <i  and  h  (defined  hy 
id.ddM  doth  exhibit  ;i  sin»;ularity.  TIkti'  is  some  correlation  with  tlu'se  results  and  those 
of  the  fully  numeric.al  schenie.  Indei'd.  it  is  most  renuirkable  th.it  the  location  of  this 
lowei  neutral  ])oint  is  identical  with  flnit  found  in  the  inconipn'ssible  casi'  ( Stewart  son  A' 
Leibovich  1987).  and  that  both  tlu' upper  and  lower  neutial  points  are  c;ipture<l  by  (0.43). 

It  is  also  worth  notinsj;  that  there  is  some  simiharity  between  the  structure  of  these 
modes,  and  those  found  in  the  incompressible  work  of  Sti'warfson  A'  Brown  (1985).  ah 
though  the  particular  det.ails  tire  different,  and  in  onr  case  aindytic/asymiitotic  solution  on 
the  /■  -  0(1)  scide  is  possible  (see  C.G).  whilst  Stewartson  A:  Brown  (1985)  ha.l  to  resort 
to  ;i  numerictd  ajiiu'oach  for  this  scale. 

In  the  followmu;  subsection  we  »;o  on  to  draw  a  number  of  conclusions  from  our  study. 


:u 


7.  Conclusions 

\V('  liav('  inouiit('(l  a  systc'iiiatir  study  of  tlu'  iiiviscid  stal)ility  of  tlu'  trailin,!;  line 
voit('x.  starting  at  Z(To  Mach  imnilx’r  M ,  and  ])rogr('ssiv('ly  incn'asing  M .  \\v  sec  a  gi'iicial 
r('duction  in  growth  rates  as  M  incn'ast's,  and  indet'd  tlu'  rt'sidts  of  suhsi'ctions  5.4  and 
5.5  predict  that  the  original  family  of  mode's  will  stahili.se'  wlu'ii  M  —  ()(|n|^  ).  specifically 
win'll  (5.43)  is  satisfic'd.  Howe've'r  it  is  shown  in  se'ction  C  that  wlu'ii  M  =  ()(|n|).  a  centre 
mode'  class  of  instahility  is  foriiK'd.  We  fe'c'l  that  although  our  numerical  if'sults  (figures 
12  in  particular)  indicate'  that  tlu'se'  mode's  spring  from  the  highe'r  orde-r  modes  at  lower 
Mach  numhe'rs.  as  —  n  incre'ase's  the'se'  ce-ntre'  mode's  mode's  may  well  he'e-ome'  elistine’t  freim 
the'  original  e-lass  of  meide's  that  e-xist  at  lowe'r  Mae-h  numhe'rs. 

A  furthe'r  impeu  tant  fe-ature'  of  neite'.  anel  enu'  that  is  eihse'rve'el  in  ine-oniiue'ssihle'  weirk 
(Le'ihewich  A'  Ste'wartsein  19S3.  Ste-wartsem  A:  Breiwn  19S5)  is  that  tlie'ie'  is  a  geioel  ele-al  e)f 
nume'rie-al  e'viele'iice'  to  sugge'st  that  the'ie'  e'xist  ne)  instahilitie's  for  o  >  —t  hi-  Inele'e-el.  this 
is  also  e-onfirme'd  hy  (G.45).  neiting  that  a  is  give'n  hy  (G.33).  This  asjie'e-f  is  e-urre'iitly  unele’r 
furthe'r  inve-stigation. 
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Figure  11c  Variation  of  growth  rate  with  M  using  section  5.4  results, 
first  6  modes. 
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Figure  12a  Variation  of  growth  rate  with  M  fully  numerical  results,  n 
first  6  modes. 


5.5  results,  n  =  —15,  a 
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Figure  14b  Variation  of  growth  rate  with  a  for  centre  modes,  n  =  —10,  M  —  10. 
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Figure  15  Variation  of  growth  rate  with  a  n  =  — 10,  M  =  10,  fully  numerical  results. 
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